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ABSTRACT 

We have searched for prompt radio emission from nine Gamma Ray Bursts (GRBs) with a 12 m 
telescope at 1.4 GHz, with a time resolution of 64 /xs to 1 s. We detected single dispersed radio pulses 
with significances > 6a in the few minutes following two GRBs. The dispersion measures of both 
pulses are well in excess of the expected Galactic values, and the implied rate is incompatible with 
known sources of single dispersed pulses. The arrival times of both pulses also coincide with breaks in 
the GRB X-ray light curves. A null trial and statistical arguments rule out random fluctuations as the 
origin of these pulses with > 95% and ~ 97% confidence, respectively, although a simple population 
argument supports a GRB origin with confidence of only 2%. We caution that we cannot rule out 
RFI as the origin of these pulses. If the single pulses are not related to the GRBs we set an upper 
limit on the flux density of radio pulses emitted between 200 to 1800 s after a GRB of 1.27?i'^^/^ Jy, 
where 6.4 x 10~^ s < w < 32 x 10~^ s is the pulse width. We set a limit of less than 760 Jy for long 
timescale (> 1 s) variations. These limits are some of the most constraining at high time resolution 
and GHz frequencies in the early stages of the GRB phenomenon. 
Keywords: gamma-ray burst: general 



1. INTRODUCTION 
When p rocessing archival dat a from the Parkes 64 m 



at lower frequencies ( ~ 100 MHz) both for for the col- 



telescope, Lorimer et aL| ( 2007 ) detected a single 30 Jy 
burst with a spectral index ot -4 {S (x v") with ~ 5 ms 
duration at a dispersion measure (DM)of 375 pc cm~^. 
The dispersion measure was much larger than the Galac- 
tic contribution in the direction ( Cordes & Lazio 2002 1 
implying an ex t ragala ctic origin. 1^'rom this detection, 
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ture around ^ 
lO-^hr-ideg- 



(^2007) 
10^ 



derived 
K and 



a brightness tempera- 
n event rate of 3.8 x 



(based on a sample of 1). Lorimer et al 



(| 2007[ ) proposed that these parameters were broadly 
compatible with a Gamma Ray Burst (GRB) origin, but 
noted that no mechanisms had been discussed that could 



produce such a burst. Keane et al. (2011) also report a 
detection with the same telescope of a somewhat weaker 
burst (4 Jy) with higher DM (745 pc cm~^). 

To date, a GRB origin for these short-timescale, GHz 
radio pulses has yet to be observationally tested. Fur- 
thermore, some doubt has been cast on the as tronomical 
jin of these bursts (Burke-Spolaor et al. 2011a Kocz 



and the state ol published theoretical mech- 
anisms with the observed timescales and frequencies has 
not progressed. There are some, as yet unobserved mech- 
anisms that can produce longer timescale radio emission 
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and short GRBs (|Lipunov & Panchenko|1996| 


Hansen &; 


Lyutikov||2001l iPskirkov k Posinov^OlUl |yiii 


bata et al. 



201 1[ ), and scattering effects may limit the observ ability 



of short-timescale bursts in certain circumstances ( Mac- 
[quart.2007 ; ,Lyubarsky 2008) . 

INonetheiess, the rewards for detecting prompt emis- 
sion from GRBs are great. Many aspects of the ex- 
plosion physics can be probed by measurements of the 
prompt radio emissi on, such as the jet opening angle, 
and Lorentz factor (|Macquart||2007) , the density an d 
distance of any scattering materiar^(Lyubarsky 2008 L 
and the str ucture of the fireball magnetic field ( Sagiv 



et al. 20041. Additionally, a short radio pulse from an 
extragalactic source is expected to undergo dispersion as 
it propagates through the intergalactic medium (IGM). 
Such dispersion would not only provide direct evidence 
for the e xistence of the m ajority of the baryons in the 
Universe ( Ginzburg||1973 ), but, for bursts of sufficiently 
high redshift, would also differentiate between diff erent 
models of cosmic reionization history ( Inoue|[2b04 ) . 

There have been a number of unsuccessful searcnes for 
prompt emission from G RBs, althoug h no searches were 
sensitive to the burst of j Lorimer^t aL| ( 2007 ) . Most 
searches have been non-directed, in which a large frac- 
tion of the sky was monitored, with the hope of a GRB 
occurring somewhere in this region. Early results at 151 
and 408 MHz with an integration time of 300 ms de- 
tected some pulses within ±10min of the gamma-ray 
trigger, but did not confidently associate a ny with GRBs 
( Cortiglioni et al!]|1981 Inzani et al.||1982 ). A survey at 
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843 MHz, which was sensitive to pulses between 0.0 01 ms 
and 800 m s, also made no definitive detections (Amy 
et al.||1989 1 , although it was not clear if any GRBs were 
present m the field of view d uring the observations. More 
recently, Katz et al. (2003) performed an all-sky survey 
at 611 Mhlz, with a time resolution of 125 ms and a flux 
density detection threshold of 27 kJy. This search de- 
tected ~ 4 X 10^ bursts in 18 months, but rejected 99.9% 
as RFI and identified the remaining bursts w ith solar ac- 
tivity. In a similar vein, Lazio et al. (2010) detected no 
transients above 500 Jy for pulse widths ot about 300 s at 
73.8 MHz. A number of other surveys have begun but are 
yet to pubhsh results ( Balsano et al.|1998 Morales et al. 
20051. The only report ot autom atic toilow-up a t radio 
treq uencies is that perform ed by Koranyi et al. ( 1995 1 
and Dessenne et al. (1996) at 151 MHz with a time res- 
olution ot 1.5 s. Based on observations of two GRBs, 



Dessenne et al. ( 1996) report upper limits on any radio 
emission of 16-73 Jy between 5 hrs before, and 2 hrs 
after the GRB. A search for the evaporation of primor- 
dial black holes at 3 GHz with a time resolution of 2/xs 
also f ailed to detect any radio emission ( O'Sullivan et al. 
19781. 



The lack of radio detections of previous surveys may 
be due to their low operating frequencies, low time reso- 
lution, insufficient sensitivity, and low sensitivity to the 
GRB rate. While low frequency observations have the 
advantage of a large predicted radio luminosity due to 
steep spectral index, which is predicted by some mo dels 
of prompt radio emission (e.g. dSagiv fc Waxman|2002|) , 
and has been observed in one case ([Lorimer et ai.||200Y' ), 
some low frequency effects make detecting short dura- 
tion pulses more difficult. For example, scatter broaden- 
ing, which subs tantially redu ces the detectability of radio 
pulses ( Cordcs fc McLaughlin 2003) , and sky tempera- 
ture are worse at low trequencies. Low time resolution 
reduces the detectability of short duration radio pulses, 
including those required to avoid the brig htness tem 
perat ure constraint from induced scattering ( L yubarsky 
20081). In p artic ular, the directed sea rches ot |Koranyi 



et al. ( 1995 ) and |Dessenne et al. ( 1996 1 used a relatively 
low time resolution ot 1.5 s. The sensitivities of most 
of the above surveys are low, and not a pproaching the 
flux le vels required to detect the burst of Lorimer et al. 
(2007), while more sensitive blind experiments have not 
had tne field of view and on-sky time to obtain a GRB 



2009) 



in-beam ([W ayth et al.|2012[ [Burke-Spolaor et al.|2011b 
Deneva et al 



In this paper we describe a survey to detect prompt 
radio emission from GRBs at 1.4 GHz, to test whether 



the bursts r eported by Lorimer et al. (2007) and Keane 
et al. (2011 ) have a GRB origin. We used a single 12 m 
dish that slews automatically to the GRB coordinates, 
based on the gamma-ray position, and observes the po- 
sition with high time resolution. Our aim was twofold: 
to attempt to detect any GHz radio emission within the 
first few minutes of the GRB, and to gain experience in 
automating radio follow-up for a potential future exper- 
iment. 



2. OBSERVATIONS 
2.1. Telescope Configuration 



We used a single 12 m dish antenna at the Parkes ra- 
dio observatory. The visible declinations ranged between 
-90deg < i5 < -h20deg. The dish had a 1.4 GHz horn at 
the prime focus, with orthogonal linear polarizations sup- 
plying room temperature low noise amplifiers bolted onto 
the horn, and wrapped in a thermally insulating blanket. 
The Tgys was measured at approximately 100 K and the 
bandwidth was 220 MHz between 1230 to 1450 MHz. 
We estimate the system equivalent fiux density to be 
3.8 X 10"^ Jy. No radiometer switching was used. Both 
polarizations were upper-sideband mixed down to a cen- 
tre frequency of 600 MHz and digitized and channelized 
by a digital spectr ometer based on an IBOB (Internet 
Break-out Board) (McMahon 2008). Sampling was per- 
formed at 800 MHz, resulting m sampling at the second 
Nyquist zone. The sampled signal of each polarization 
was channelized with a 1024 channel polyphase filterbank 
spread across 400 MHz, with an approximate frequency 
resolution of 390 kHz. There were approximately 560 us- 
able channels across the 220 MHz bandwidth. The cross 
polarization product was not formed. The channelized 
output was detected, integrated and dumped every 64 
microseconds with a resolution of 8 bits per channel. The 
integrated spectrometer data were written to disk in real 
time and all processing was performed offfine. 

2.2. GRB response 

A control computer was connected to the Gamma-ray 
Coordin ate Network (GCN) via the socket distribution 
method ( [Barthelmy et al.||1995| . The GCN was config- 
ured to send only notifications when the time delay be- 
tween the satellite detection and socket distribution was 
less than 1 hour, the position of the burst was more than 
10 degrees above the visible horizon, and the position er- 
ror was less than 1 degree (the width of the primary beam 
of the telescope). The GCN was configured to send noti- 
fications from the Swift /Ji AT, Fermi, Integral and Agile 
missions, although only events from Swift /BAT passed 
the filter constraints. 

When a GRB alert arrived from the GCN, the tele- 
scope immediately began saving the 8-bit channelized 
data. It then slewed to the coordinates given in the burst 
alert, while saving the raw spectra to disk. Typically the 
slew took two minutes, after which the telescope was "on 
source" . After 30 minutes on source on the initial GRB 
position (the tracking position was not updated if addi- 
tional position refinements were sent by the GCN), data 
capture was stopped. To check that the receiver sys- 
tem was operating correctly, immediately after the GRB 
capture stopped, the t elescope slewed to a bright pulsar, 
eith er the Vela pulsar ( iLarge et al.|1968| ) or PSR B1641- 



45 (Komesaroff et al. 



1973| and took an additional 10 



minutes of raw data once the antenna was on source 

The telescope performed pulsar monitoring when it 
was not observing a GRB alert, or observing a bright 
pulsar for the system check. The pulsar monitoring re- 
sults are not reported here. 

2.3. Processing 

Our primary aim was to search for repeating signals, 
and single pulses at high time resolution (~ 1 ms). Our 
secondary aim was to search for slowly varying signals 
at low-time resolution (> 1 s). For the slow search, dis- 
persion can be ignored as the dispersion delay across the 
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observing bandwidth is ^ 770 ms for a dispersion mea- 
sure (DM) of 1000 pc cm~'^, which is less than the time 
resolution. Both timescales were affected by RFI, which 
can increase the number of false detections. 

2.3.1. Kurtosis-based RFI excision m the frequency domain 

When a channel is affected by impulsive RFI, it typ- 
ically contains more samples with large values than a 
channel containing Gaussian noise alone. One method 
for measuring an excess of large values is to compute 
the 'excess kurtosis' (k), which compares the number of 
large values in a time series against the expected number 
of large values given a Gaussian distribution. 

To determine which of the channels contained impul- 
sive RFI, we took the 8-bit raw data and, for each channel 
computed excess kurtosis in one second intervals. If we 



define 



as the n sample in a given interval, and N 



as the number of samples in the interval, then the excess 
kurtosis is given by the equation: 



„ 



(1) 



where ^4 is the fourth moment of the distribution, de- 
fined as 



JV-l 



M4 = 



n=0 
N-1 



-E 



n=0 



(2) 
(3) 



and cr"* is the square of the variance of the distribution, 
given by 



A^- 1 



E 

n=0 



{x{n\ - iif 



(4) 



The excess kurtosis is calculated for each channel in 
one-second bins, and displayed in an interactive program. 
This program plots the maximum and average kurtosis 
as a function of channel number. The maximum kur- 
tosis for a channel that was not affected by RFI was 
K ~ 0.6, and the RFI affected channels had larger values 
of K (Figure [ij . An RFI channel mask was created by 
setting all channels with an average k, less than a thresh- 
old to a weight of 1, and the remaining channels to a 
weight of zero. The threshold could be set interactively 
in the program, and was set to between 0.6 and 0.8. In 
a subsequent step, a user-defined number of additional 
"padding" channels were set to zero on either side of a 
channel that had zero weight because of the threshold 
step. The number of additional padding channels was 
between 1 and 3. 

The kurtosis was calculated separately for each polar- 
ization, but the same threshold and number of padding 
channels were used for each. This resulted in different 
channel masks for each polarization. The final channel 
mask was the logical 'and' of both channel masks. 

The resulting channel masked was used for both for 
low-, and high-time resolution processing. 



2.3.2. Low-time resolution processing 

For low-time resolution processing, we took the mean 
of each channel in one-second bins. To compute a light 
curve, we weighted the spectra by the channel mask that 



we calculated according to the method in section 2.3.1 
and took the mean value across frequency. We then took 
the mean of the light curves of both polarizations. 

The final product was a radio light curve, which we 
visually inspected for variability or transient emission. 

2.3.3. High-time resolution processing 

Each channel in each polarization was offset to have 
zero mean and scaled to a variance of unity. The offset 
and scale were updated every 10 seconds. Corresponding 
channels in both polarizations were then added together 
and the result truncated to 2 bits. A periodicity search 
was perform ed in an almost identical manner to |Keith| 
et al. (2010 ) and single- pulse search was performed iden- 
tically to Burke-Spolaor et al. (2011b). Each time se- 
ries was dedispersed with 1991 DM trials between and 
lO'' pc cm~'^. The maximum DM was chosen to be well 
in excess of the expected Galactic and host contribu- 
tions (probably < 1000 pccm"'^ in both cases), and the 
predicted intergalactic medium for GRBs out to redshift 
z ~ 6 of > 6000 pc cm-3 (loka 2 003] |lnoue||2004| ). 

DM steps were chosen so that the smearing due to the 
finite DM step was at most 1.25 times the DM smear- 
ing in a frequency channel at the center frequency. No 
resampling for DMs larger than the diagonal DM was 
performed. This strategy results increasing DM steps, 
with small spacings at lower DMs and larger spacings 
at higher DMs. Time- based RFI excision was applied. 
Unlike the processing of Keit h et al. (20101, a frequency 
mask (calculated in Section 2.3.1 ) was applied, instead 
of the fast Fourier transform (FFT) based channel mask- 
ing approach. An FFT based harmonic search was per- 
formed to search for periodic or pulsar-like signals, and 
a single-pulse search was also performed to search for 
bursty signals. The single-pulse search used the dedis- 
persed time series, and applied 9 different boxcars of 
width 2* samples (0 < i < 9) to search for pulse widths 
ranging from 64^s to 32 ms. The single-pulse candidates 
were grouped across DM and boxcar trials using a friend- 
of- friends algorithm (Burke-Spolaor et al. 2011b). This 
algorithm exploits the tact that a broad- band pulse well 
above the minimum detection threshold appears in mul- 
tiple DM and boxcar trials. Statistical variations near 
the threshold tend to appear in only a single DM and 
boxcar trial. Therefore, this algorithm declares a candi- 
date only if a pulse is detected in a group of 3 or more 
adjacent DM and boxcar trials with each detection hav- 
ing a significance of > 6a. Groups with less than three 
detections were classified as Gaussian fluctuations, and 
groups whose peak significance was at a DM of less than 
2 pc cm~^ were classified as RFI. 

Visual inspection of the resulting candidates clearly 
showed that our sample contained substantial impulsive 
RFI. In particular, we observed impulsive RFI that was 
detected strongly at low DMs and then disappeared at 
intermediate DMs, only to reappear again at higher DMs. 
The time difference between high-DM and low-DM 'is- 
lands' of detections was occasionally as large as 2 s. The 
exact cause for this disappearance is not clear. The result 
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was that the low-DM detections were classified as RFI 
(as they peaked near zero DM) , but nearby high-DM de- 
tections were erroneously classified as detections, as the 
friend-of-friends algorithm failed to join a continuous- 
stream of detections from the high DMs to the low DMs. 
To counteract the island problem, we rejected high-DM 
candidates that had time-coincident low-DM detection. 
I.e. we re-classified candidates initially produced by the 
pipeline as RFI if there was a > 6(t detection within 3 s 
of the candidate and with a DM of < 25 pc cm~^. Visual 
inspection of the resulting candidates confirmed that all 
remaining RFI was excised by this additional filter, and 
that no believable low-DM candidates were excised. 

Finally, we visually inspected plots of the candi- 
dates which were pro duced by the pipe line using the 
PSRCHIVE package ( |Hotan et ar]|2004[ ) and our own 
custom tools. 

3. RESULTS 

We describe the results of the single pulse, periodic 
and low-time resolution searches in this section. We note 
that the dispersion delay between the gamma-rays and 
the 1.4 GHz observing frequency, at the maximum DM of 
10000 pc cm^-^ is < 30 s, which is small enough to ignore 
in the following analysis. 

3.1. Events 

The telescope was operating between 2010 June 
and 2011 February, during which it responded to 16 
Swift /BAT alerts of which nine were confirmed GRBs 
(Tablejl]). The remainder were instrumental effects, asso- 
ciated with known recurring objects, or Galactic sources 
such as X-ray binaries. 

Of the nine GRBs, two had a duration < 2 s, satisfy- 
ing the definition of short GRBs. The remaining seven 
GRBs were long GRBs. We detected no radio counter- 
parts in the single pulse, periodic or long integrations for 
any short GRBs. 

3.2. Single-pulse search 

The friends-of-friends algorithm detected candidates 
for all GRBs we observed, and after the post-filtering 
step was applied, eleven candidates remained. We re- 
jected nine of these candidates because they had rela- 
tively low DM (< 200 pccm""^) and occurred during pe- 
riods of frequent, strong RFI (i.e. within ^ 100 s, of the 
candidate). 

Of the nine GRBs for which we obtained data, eight 
had reasonably interference-free observations (the excep- 
tion being GRB 101020A, which we will not consider 
any further, due to substantial flagging; see Fig. [l] , 7th 
panel) and two GRBs had one single radio pulse which 
was not associated with low-DM RFI (Figure [2]) . For 
GRB 100704A, a 6 ms wide single pulse was detected 
1076 seconds after the GRB at a DM of 195 pc cm"^, at 
a significance of 6.2a . This pulse is fairly indistinct in 
the time domain (Figure[3|. For GRB lOlOllA, a 25 ms 
single pulse was detected 524 seconds after the trigger at 
a DM of 570 pccm"'^, at a significance of 6.6cr (Figure 

i- 

The low significance of both detection means that fit- 
ting a spectral index is difficult. Both pulses appear to 
be fairly consistent with a spectral index of zero, which 
we will assume for the remainder of the analysis. 



Both pulses have approximately the same energy in 
both polarizations (Figure [sj |4| , arguing that they are 
either unpolarized, circularly polarized, or linearly polar- 
ized with a position angle approximately 45° to the feed 
angles. A linearly polarized pulse parallel to the feed 
axis can be definitely ruled out by our data, as such a 
pulse would appear in only one polarization. The chance 
of the polarization angle being aligned at 45° to the feed 
is small. We therefore consider it more likely that the 
pulse is unpolarized, or circularly polarized. 

3.2.1. Event rate 

We detected two pulses in 4 hrs of observing. The full 
width, half maximum beam width of the 12 m antenna 
at 1.4 GHz is approximately 1.3 deg, which implies an 
area of 1.3 deg^. If these events are not related to the 
GRBs and could have been discovered in a blind search, 
then the implied event rate is 2.9 x 10~^ deg^^hr^^. 

3.2.2. Correspondence with X-ray light curves 

Swift/XRT light curves are available for both GRBs 
for which we detected si ngle pulses. For GRB lOlOllA, 
Cannizzo et al. (2010a I fit a three component broken 
power law to the 0.3-10 keV light curve, and derive a 
value of t he final break in the po wer law of tbrcak,2 = 

do not describe their 
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+ 175 , 
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Cannizzo et al. 



(2010a 



exact fitting method. We nt 



the same data with a broken 



power law, with each piece having a form S(t) cx t"^ for 
< i < iVbreak, and with each data point weighted by 
the flux error. We set the ini tial break times a.t t = 
116 and t = 602, as derived by Cannizzo et al. (2010a). 
From this fit (x^/A'd.o.f = 1-6) we derive a break time of 
ibrcak,2 = 707 ± 173. The arrival time of the radio pulse 
coincides with the break in the power law, to within the 
statistical errors for both fitting methods (Figure [5|. We 
note that the dispersion delay between the gamma-rays 
and our radio frequencies is less than 3 s and therefore 
insignificant for this calculation. 

For GRB 100704A, no fi t is derived in the GCN re- 
port (Grupe et al. 2010b I. We fit a four component 
power law to the A-ray light curve for the data after 
t = 400 s (Figure[5]), corresponding with the beginning of 
the proportional-counting mode data set. We used initial 
breaks times at 500,1000 and 3 x 10^ s. The resulting fit 
ix^/Nd.o.f = 1-2), has a break at t = 1700 ±410 s, which 
is within 1.5a of the pulse arrival time a.t t = 1076 s. We 
note that the time of the fitted breaks is quite sensitive 
to the choice of initial break times. 

3.3. Periodicity search 

Our method detected candidate repeating signals for 
all GRBs we observed, however, visual inspection of the 
candidate plots indicated that almost all of these candi- 
dates were associated with RFI. 

The pulsar search yielded a single candidate for a re- 
peating signal that was not conclusively associated with 
RFI. The candidate was detected in the data following 
GRB 110412A with a S/N of 5.3 in the periodicity search, 
a topocentric period 524.48 ms and a DM of 51 pc cm~^. 
After optimization for DM and period, the S/N of the 
candidate increases to 6.1. The candidate is brightest 
during the first 20 minutes of the observation. The sig- 
nal is clearly detected in all sub-bands, as expected for 
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Figure 1. Plots of time- aver aged long-time integrations for each GRB observation. For each GRB top panel: bandpass, 2nd panel: 
standard deviation, 3rd panel: kurtosis, 4th panel: channel weights. Solid lines are the mean values over 30 min, and dotted lines are the 
maximum values over 30 min. 
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Figure [T] continued. 

a broad-band pulsar signal, and which also rules out 
narrow-band RFI. The profile comprises a single com- 
ponent with a width of about 5 ms . The signal is not 
very well defined in the period-DM plane or in the DM- 
S/N plot, but this is not unusual as the S/N is relatively 
low. 

This position has been surve yed by the Parkes 70 cm 
survey ( Manchester et al.|1996 ) with no detection. Given 
our sensitivity, any pulsar detected by our experiment 
would have been detected in that survey with very high 
significance, u nless the pulsar was nulling at the time 
( Backer|1970 l, or the detection was a transient signal re- 
lated to the GRB onset. It is common for low S/N candi- 
dates which appear convincing the first detection to not 
subsequently be detected in follow-up (M. Keith, private 
communication). On these grounds, the astronomical na- 
ture of this detection is very difficult to confirm and we 
will therefore not consider this candidate any further. 



3.4. Low-time resolution search 

The light curves at 1 s integration are shown in Figure 
[6] The light curves with rapid changes in flux density 
can be attributed to interference. For example, the flux 
density variations during observations of GRB 1000823A 
are due to RFI from the Beidou Gl satellite|£] This 
can be seen from the modulation spectrum centered at 
1270 MHz in the bandpass and variance plots. During 
the observations, the telescope pointing direction was 
within 10 degrees of the satellite position when the in- 
terference was most severe (Figure [7f. Ahhough our 
method flagged the worst channels, it is clear that ad- 
ditional power leaked through the unflagged channels, 
affecting the light curve. 

The substantial flux variations during observations of 
GRB101020A can be attributed to the proximity of the 
Sun to the observing direction, which was 25 degrees 
from the Sun at the beginning of the observation. None 

^ Refereed, English information about the Beidou global satel- 
lite system is non-existent. The only available information avail- 
able in English is on general news and satellite navigation web- 
sites. Perhaps the best source of information can found at http: 
|77en.wikipedia. org/wiki/Beidou_navigatlon_system The otti- 
cial Uhinese website is http://www.beidou.gov.cn Beidou satel- 
lite ephemerides are provided by the US Space Command and the 
real time position was computed using data provided at ,http : | 
/ /www . n2yo . com/ 



of the narrow-band RFI coincides with the large changes 
in flux density, indicating that the changes in flux are 
broad-band and are likely to be a result of the Sun mov- 
ing through the far sidelobes. 

The variations during observations of GRB 110412A 
are likely due to unflagged narrow-band interference at 
around 1230 MHz (Figure |8]). These channels had no 
substantial kurtosis and so were not flagged, but had 
slightly higher variance than the rest of the band (Figure 
fl]). Nonetheless, the fact that the interference is narrow 
oand strongly suggests RFI as the source of the variation. 

The slow variations (e.g. GRB 100625A) are due to 
gain variations in the system, which is stable to about 5% 
in 30 minutes. We were unable to calibrate out these gain 
variations, as our feed was not equipped with a switched 
radiometer. 

4. DISCUSSION 

4.1. Single-pulse detections 

We have two single-pulse candidates, the brightest of 
which is a 6.6ct pulse 524 s after GRB lOlOllA. Below 
we discuss possible origins for this pulse. Our conclu- 
sions also broadly apply to the other candidate single- 
pulse (the 6.2cr pulse 1076 s after GRB 100704A), but 
we restrict our attention to the brightest pulse. 

The DMs of both pulses are much larger than the 
DM in the GRB direction pr edicted by the NE2001 



|2002[ ) for our Galaxy. For 
NF20U1 predicts a maximum DM of 



model (Cordes & Lazio 
GRB lOlOUA, 

39 pccm"'^ at a DM distance of 3.14 kpc, which im- 
plies an excess DM of 530 pccm^^. For GRB lOlOllA, 
NE2001 predicts a maximum DM of 106 pccm-3 at a 
DM distance of 9 kpc, which implies an excess DM of 
89 pccm^'^. The large excess DMs for both radio pulses 
are well above the uncertainties in the NE2001 model 
and therefore, if the pulses are indeed astronomical, ex- 
clude Galactic sources such as pulsars or RRATs. The 
excess DM must therefore be due either to intergalactic 
material or a possible host galaxy. 

4.1.1. Random noise fluctuation - theoretical arguments 

In principle we can apply statistical arguments to de- 
rive the probability of a 6.6cr pulse occurring by chance in 
our data. In practice, such arguments are not straight- 
forward for a number of reasons. Firstly, our data are 
not easily described by analytic probability distribution 
functions, chiefly because of the presence of RFI. Sec- 
ondly, the processing pipeline applies a number of RFI 
mitigation strategies and matching algorithms that are 
difficult to model in a statistical sense. Finally, the search 
over a range of DMs and pulse widths, which necessar- 
ily sample overlapping regions of parameter space, is not 
independent. That is, a genuine astronomical signal will 
appear in a number of adjacent DMs and boxcar trials, 
implying that the simplifying assumption of independent 
trials is not valid. 

Nonetheless, we can construct a simple model with 
which to evaluate our real-world results, assuming Gaus- 
sian noi se and independent trials. Here we follow the 



logic of Burke-Spolaor et al. (201 lb I. For a 30 minute 



observation, we produce 1991 DM trials, each compris- 
ing Ns = 28.125 X 10^ samples. We apply a set of 9 
boxcar trials to each DM trial, implying that a total of 
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Figure 2. DM vs. time for the two GRBs with single-pulse candidates. Single-pulse detections with a significance > 6(t appear as circles 
in this plot, with the size proportional to the S/N. The detections are color coded according to their classification by the friends-of-friends 
algorithm as candidates (blue), false positives (red) and RFI (grey). Top panel: GRB 100704A with a single pulse candidate 1076 seconds 
after the GRB at a DM of 195 pc cm~^, with a significance of 6.2a and width of 6 ms. Bottom panel: GRB lOlOllA with a single pulse 
candidate 524 seconds after the GRB at a DM of 570 pc cm"'^, with a significance of 6.6(t and width of 25 ms. The time origin of these 
plots is the time that the telescope first arrived on source (Ton). For clarity, DMs from 1000 pc cm""^ are not shown. 

These relatively high probabilities suggest that it is 
quite likely that we would find a 6.6a candidate in a 
single observation, and that we would also expect several 
candidates across our 9 observations, exactly matching 
our results. However, we have neglected up to this point 
the 'friend-of-friends' algorithm employed by the single- 
pulse pipeline, which declares a candidate as viable only 
if a pulse is detected in 3 or more nea rby DM trials. 
According to Burke-Spolaor et al. (2011b), this algorithm 
reduces the number of candidates by a factor of ~ 100 — 
1000, implying that applying the algorithm reduces the 
probability of detecting exactly one candidate from 4.6% 
to 4.4%. Such an analysis suggests that the 6.6(7 pulse 
we have detected is not a random noise fluctuation, with 
a confidence of greater than 95%. 

Relaxing the assumption of independent trials implies 
that adjacent trials (either in DM of boxcar width) have 



Np = 1991 X ^^^0 N,/2' ~ 1.1 X 10" points are searched 
for pulses. Assuming each trial is independent, then the 
expected number of false alarms exceeding a threshold T 
during an observation is A^perfc(r/(cr\/2)), where erfc is 
the complementary error function. 

Applying this model to observations of GRB 101011 A, 
we expect Nqo- — 216 trials above 6a and get 260 such 
detections, roughly tallying with our discussion to this 
point, but many detections are clearly associated with 
RFI and are not related to the Gaussian distribution. At 
the 6.6(7, significance of our detected pulse, we expect 
on average 4.6 false alarms per observation. Applying 
the Poisson distribution, we can compute the probability 
of finding exactly one false alarm during a 30 minute 
observation, at or above 6.6a, of 4.6%. Alternatively, 
we can compute the probability of finding one or more 
pulses above 6.6a of 94%. 
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Figure 3. Detail of tiie 6.2(t pulse detected 1076 seconds after 
GRB 100704A. The top three panels are the dedispersed time series 
of 4 frequency channels spread across the band. The pulse appears 
clearly in both polarizations (the top 2 panels) and in the sum of 
the two polarizations (third panel). The bottom panel is the time 
series where all the frequency channels have been summed. The 
origin of the time axis is the pulse arrival time (Tpuigc). 



overlapping time and frequency samples. Thus, if a given 
trial has a high value, an adjacent trial is also likely to 
have a high value, because many of the samples overlap. 
Therefore, relaxing the assumption of independence in- 
creases the false alarm probability, and we would expect 
to see clusters of false alarms in DM and boxcar space. 
The amount of overlapping samples between adjacent tri- 
als is of the order of 50% both in time and frequency. For 
a trial with a 6<t detection, which is already unlikely, the 
probability of the non-overlapping samples from an adja- 
cent trial pushing the adjacent trial above 6a is small. As 
the friend-of-friends algorithm requires three or more ad- 
jacent trials above 6a to declare a detection, it effectively 
nullifies the correlation between adjacent trials, and the 
probabilities described above are still accurate. 

For simplicity, we have neglected the effect of the ad- 
ditional weak RFI filter, which discards can didate s with 
adjacent low DM pulses within 3 s (Section [2.3.3 ). Such 
a filter will reduce the number of candidates produced by 
random fluctuations even more, which in turn increases 
the confidence. 
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Figure 4. Detail of the 6.60- pulse detected 524 seconds after GRB 
lOlOllA. The top three panels are the dedispersed time series of 
4 frequency channels spread across the band. The pulse is clearly 
detected across all channels, in both polarizations (the top 2 pan- 
els) and in the sum of the two polarizations (third panel). The 
bottom panel is the time series where all the frequency channels 
have been summed. The origin of the time axis is the pulse arrival 
time (Tpuise). 

4.1.2. Random noise fluctuation - A null trial 

The theoretical considerations in Section |4. 1 . II assume 
Gaussian noise, and do not deal with the complexities of 
the RFI and detection algorithms very adequately. Per- 
haps a more reliable approach to determining the prob- 
ability of such a pulse occurring by chance is to measure 
the number of false positives recorded when there is no 
GRB in the beam. One such method would be to ob- 
serve a 'blank' patch of sky, not containing a GRB, for a 
period of time, and measure the number of candidates. 
Unfortunately, such observations were not made when 
the antenna and feed were available. 

Ideally, we would like to remove the astronomical sig- 
nals from the data, while preserving the noise proper- 
ties and RFI, so that an accurate measurement of the 
background can be performed. Once such approach is to 
search for pulses in the data we already have, but ran- 
domize the channel ordering before dedispersion. Ran- 
domizing the channels destroys the t cx dispersion 
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Figure 5. The XRT (0.3-10 keV) light curves of the two GRBs 
after which we detected a single radio pulse. The arrival time of 
the radio pulse (Tpuige; see Tablelll is marked with a vertical line. 
Top panel: the X-ray ligh t curve of GRB lOlOllA . The black line 
is fit is the fit derived by |Canmzzo et al.| (|2010a[ l. Bottom panel: 
the X-ray light curve of (JKB iUYU4A. ihe black line is the fit 
we derived. WTSLEW is the windowed timing mode during the 
slew of the spacecraft, during which flux measurements can be less 
reliable. WT is windowed timing mode. PC is photon counting 
mode. 

relationship expected from an astronomical pulse, and 
almost guarantees that no genuinely astronomically dis- 
persed pulses will be visible in the results, especially at 
large DMs. Randomizing the channels still preserves the 
noise properties, and narrow-band, and broad-band im- 
pulsive RFI, as desired, but impulsive RFI with interme- 
diate bandwidth is destroyed. The level of intermediate 
bandwidth RFI is therefore reduced in the randomized 
data, and the number of background detections is con- 
sequently reduced. We expect the amount of intermedi- 
ate bandwidth RFI to be small in comparison with the 
broad-band RFI, so we expect this effect to be small. 

This technique, while inferior to taking additional data 
on a blank patch of sky, still has advantages over the 
theoretical approach. In particular, it uses the same 
detection algorithms as those used on the on-sky data, 
the noise statistics of the data are preserved, and the 
zero DM sequence (which is used for RFI flagging by the 
pipeline) is also preserved. 

We conducted such a trial based on the data from GRB 



101011 A (which contained the brightest single pulse). 
A software-based pseudo-random number generator was 
seeded with random data from the Linux kernel entropy 
pool at the beginning of each observation, and the or- 
dering of the channels was changed randomly every 25 s 
during the observation. 

We simulated 248 30 minute observations with the ran- 
domized channel ordering, which yielded seven single- 
pulse candidates. From this trial, we can conclude that 
the probability of detecting a > Gcr candidate that passes 
all the additional criteria is approximately 3%, and we 
therefore rule out the possibility of our pulse being due 
to random fluctuations with a confidence of 97%. 

4.1.3. Impulsive RFI 

Impulsive RFI has a number of attributes that make 
it often distinguishable from astronomical sources: it is 
usually most visible at low DMs, although it can extend 
to higher ones, depending on the exact waveform shape. 
Impulsive RFI can often be clustered in time, with many 
RFI events occurring frequently over periods of minutes, 
often followed by similar periods of no RFI. For example, 
the observations of GRB 101011 A have a single RFI event 
in the first 15 minutes observing, and 38 RFI events in 
the second 15 minutes (Fig. [2|. 

The pulse seen just after uRB lOlOllA has a very 
high DM of 570 pc cm^-^, which argues against an RFI 
origin. The nearest low DM detection (a > 6(7 isolated 
pulse) appears to be some 120 s before the candidate 
pulse, which argues against the candidate being related 
to a generally high level of RFI at the time, or being a 
high-DM component of RFI which is also detected at low 
DMs (Figurej2]). 

The high DM, absence of contemporaneous low-DM 
RFI add credence to a non-RFI origin. However, in the 
absence of a coincidence detection at a geographically 
distant site, we cannot definitively rule out RFI as the 
origin. 

4.1.4. Larimer hurst events 

Our p ulse has similar prope rties to the Lorimer burst 
(LB) of 
higher t 



Lorimer et al. 



( |2007[ ). The DM of our pulse is 
;han the LB, but this can be explained by differ- 
ences in Galactic electron content, redshift and ionized 
material near the source. The Lorimer burst is shorter 
than our pulse (5 ms for the LB as opposed to 25 ms for 
our pulse). This shorter duration cannot be explained by 
reduced scatter broadening from a lower electron column, 
implying that the difference in pulse duration is intrinsic. 
The flux densities are roughly comparable (30 Jy for the 
LB, 7 Jy for ours). If the LB is associated with a GRB 
then it is likely that we have observed the same physical 
mechanism. 



In 



4.1.5. "Perytons" 
an archival search of Parkes Multibeam archive. 



Burke-Spolaor et al. (2011a I discovered 16 single pulses 
with similar spectral characteristics to those expected 
from a dispersed extragalac tic pulse. These pulses, that 
Burke-Spolaor et al. ( 2011a[ ) call 'perytons', have a num- 
ber of unique properties. In general, the perytons fol- 
low a, t (X v~'^ frequency sweep, although some pulses 
have significant departures from this sweep. The sweep 
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Figure 6. Low-time resolution 1.4 GHz light curves for all GRBs. The vertical dashed line indicates when the antenna was on source, 
with the interval preceding this line being measurements during the slew. 
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Figure 7. The raw receiver bandpass during observations of GRB 
100823A, as a function of time and frequency. The gray scale is an 
arbitrary Unear scale of the intensity. Long time integration mea- 
surements were badly affected by RFI from the Beidou Gl satellite, 
despite zero weights being applied to the obvious interference. The 
interference was centered on 1270 MHz, the transmitting frequency 
of the satellite, and the interference was worst when the pointing 
direction passed closest to the satellite. 



Figure 8. The baseline-subtracted receiver bandpass during ob- 
servations of GRB 110412A, as a function of time and frequency. 
Long time integration measurements of were affected by narrow 
band RFI at around 1230 MHz. The narrow-band interference was 
worst at around t = 1300 s, which corresponds to a bump in the 
broadband lightcurve. 

times across the 288 MHz band centered at 1374 MHz 
are between 200 and 400 ms, but mostly clustered in the 
range 320-400 ms. The perytons are extremely bright 
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(between 0.1 — 272 kJy), and have substantial amplitude 
modulations. All perytons were detected during the day, 
in general mid-morning, and mostly in the mid-winter 
months of June and J uly. The peryton event rate is 



2.3 X lO^'^deg-^hr-^ iKocz et al.l (12012') find that the 



perytons are most likely local tib'L, as they tend to arrive 
in the second half of a human second, and arrive with a 
characteristic time gap of ~ 22 s. 

The pulses we detect have somewhat different prop- 
erties. Unfortunately, we do not have sufficient S/N to 
detect small departures from the cold plasma dispersion 
law, nor any substantial amplitude modulations. Fig- 
ure |4] supports a standard cold plasma frequency sweep, 
and constant amplitude, although the significance is not 
very high. The flux densities of our pulses are approxi- 
mately 7 Jy, substantially fainter than the perytons. If 
we extrapolate the swe ep time of our br ightest burst for 
the configuration of Burke-Spolaor et al. ( 2011a) between 
1230 and 1518 MHz, we obtain a sweep time of 536 ms. 
This sweep time is substantially higher than that of any 
of the known perytons, but within the general range, 
given the small number statistics. Our brightest pulse 
was detected on 2011 Oct 11, well outside the usual Pery- 
ton season, and at 22:58 (local time) in the evening, out- 
side the usual Peryton time. There had been no rainfall 
and no lightning detected within 50 km of the telescope 
within 24 hrs of the GRB time, disfavoring the postu- 
lated lightning related cause. Finally, the event rates 
are vastly different, with our events being more common 
by over six or ders of magnitude, although |Burke-Spolaor| 
et al. (2011a) note that the Peryton event rate is highly 
non-uniform, with 12 perytons detected within 0.5 hr in 
one instance. 

From this analysis we conclude that the substantially 
higher event rate, longer sweep time, arrival date, time 
of day, and prevailing weather conditions of our pulse all 
make it unlikely to be a peryton. 

4.1.6. Gamma-ray bursts 

Two out of eightn (25%) of our observations contain 
a 1.4 GHz radio pulse, which begs the question, are the 
GRBs themselves the source of these pulses? Unfortu- 
nately, there is no known physical mechanism for pro- 
ducing such a burst within minutes of a GRB, in spite of 
the substantial theoretical work inspired by the detection 
of the LB. 

Any proposed mechanism must explain a range of 
properties. First, our two detections were only present 
following long gamma-ray bursts, so a mechanism must 
fit within the confines of the coUapsar paradigm ( Mac- 
Fadyen & Woosley 1999) (see below). The pulses are 



detected with flux densities of the order of 7 Jy, which re- 
quires a very powerful mechanism if placed at the cosmo- 
logical distances of GRBs. The short pulse width, of the 
order of 25 ms, requires a relatively compact causative 
region. Our pulses are 524 and 1076 seconds after the 
gamma-ray triggers, much later than the ~ 2 s disper- 
sion delay, so any mechanism must explain this delay. Fi- 
nally, any mechanism must also explain why both pulses 
are coincident with a break in the X-ray light curve, and 
also must explain the complete lack of flaring at X-ray or 

^ The observation of GRB 101020A is not considered in these 
statistics, due to RFI. 



gamma-ray wavelengths simultaneous with the emission 
of the radio pulse. 

Unfortunately, no host galaxies were observed for ei- 
ther GRB, so no redshift is available. For the remain- 
ing analysis, we assume a redshift of z — 2.8, which 
is approximately the m ean redshift of GRBs from Swift 
( [Jakobsson et al.||2006[ ). 

We now try to interpr et our pulses within the c ontext 
of the collapsar model ( MacFadyen k Woosley| |l999), 
namely: the collapse of a massive ( > 10 M(t, ) star pr oduc- 
ing a central engine, either black hole, or magnetar (Usov| 
1992 1, of mass ~ 2 — 3 M©; which in turn accretes mate- 
rial from a disk which produces jets, the prompt gamma- 
ray emission, and the shocks that pass through the sur- 
rounding material and into the interstellar medium. 

Concentrating on the brightest of our radio pulses, 
the flux density of approximately 7 Jy implies a spec- 
tral luminosity of 1.1 x 10^^ erg s~^ Hz~^ at redshift 2.8 
{Ho = 71kms-iMpc-\ = 0.27, Ovac = 0.73) as- 
suming a spectral index of zero, or a total luminosity of 
1.1 X 10^^ erg over a bandwidth of 1 GHz. The rest- 
frame pulse width is equal to about 7 ms (neglecting 
intergalactic scattering), implying a total energy release 
in the radio burst of 7 x 10**^ erg over a bandwidth of 
1 GHz. The energy in the radio burst fits well within 
the canonical 10^^ erg total energy release of the collap- 
sar model, implying that the energy budget can support 
such a pulse. The pulse width of 7 ms implies a causative 
region of size < 2.1 x 10* cm, some 10-100 times larger 
than the innermost stable orbit of a black hole, depend- 
ing on its mass and spin, and well within the magneto- 
sphere of a magnetar. In either case, this suggests that 
the emission emerges from a region close to the central 
engine. 

Perhaps the two most confounding problems are the 
time delay between the gamma-ray trigger and the emis- 
sion of the radio pulse, and that only one pulse is ob- 
served per GRB. If this is not an observational effect due 
to our detection threshold, the single detection implies 
that the radio pulse is related to a singular event that 
occurs 500-1000 s after the gamma-rays are released. If 
the central engine is a promptly formed black hole, it 
forms just before the release of gamma-rays. Therefore, 
there is no obvious mechanism related to the black hole, 
or GRB phenomenon, that can be invoked to explain 
the del ay. On the other hand, for a magnetar central 



engine, Baumgarte et al. (2000) suggest that a super 



Chandrasekhar (> 1.4 M©) magnetar can be rotation- 
ally supported, and can avoid turning into a black hole 
as long as its angular momentum exceeds a threshold. 
If sufhcient angular moment is lost (e.g. due to winds), 
a blac k hole is formed after some delay. [Metzger et al.| 
(2011) suggest that the timescale for delayed black hole 
formation can be a few hundred seconds, but add that 
there is considerable uncertainty in this timescale. The 
exact mechanism for how this formation could release the 
radio pulse with the observed properties is still uncertain, 
but the delayed formation model has the desired charac- 
teristics: the delayed black hole formation happens only 
once, and occurs some time after the GRB trigger. 

The alignment of the radio pulse with breaks in the 
X-ray light curve could also be related to the delayed 
formation of the black hole. The X-ray light curve is typ- 



Limits on prompt, dispersed radio pulses from gamma-ray bursts 



13 



ically explained by a number of segments ( Gehrels et al 
20091 and references therein): the initial rapidly faUing 



afterglow is followed by a plateau lasting a few thousand 
seconds, which is in turn followed by a classical afterglow. 
The so-called plateau is usually explained as being due to 
additional energy being injected into the external shock, 
for example, from magnetar winds. The radio pulse af- 
ter GRB 100704A occurs at the end of the X-ray plateau, 
which implies it is related to the end of the energy in- 
jection driving the plateau. One could imagine that the 
formation of the black hole event horizon encompassing 
the magnetar and its winds could be responsible both for 
the radio pulse and for the end of the X-ray plateau. On 
the other hand, the radio pulse following GRB 101011 A 
occurs at the beginning of the X-ray plateau (which is 
not as clear as the plateau in GRB 100704A), casting 
some doubt on the delayed black hole formation as the 
source of this pulse. 

The coincidence with the X-ray breaks could also be a 
statistical fluke. For GRB lOlOllA, the formal error in 
the break time is roughly 200 s, which is 10% of the ob- 
servation duration. There are at least two breaks in the 
X-ray light curve predicted durin g the observation inter- 
val of 200-1800 s ( [Gehrels et ar|f2 009), suggesting a 20% 
chance of a random event tailing within the statistical er- 
rors of a break in the X-ray light curve. The constraints 
on the break time in GRB 100710A are less strong, im- 
plying an even higher chance of a random event falling 
within the errors of an estimated X-ray break. 

4.1.7. Not a GRB 

Given the spatial distribution of GRBs, we would ex- 
pect to receive pulses with a range of flux densities, both 
near our threshold and far above it. The key argument 
against a GRB origin for these pulses is the small like- 
lihood of our detections being so close to our detection 
threshold and the absence of any more significant detec- 
tions. 

To quantify this effect, we performed a monte-carlo 
simulation of 5 x 10'^ GRBs. Our simulation assumed 
the GRBs were distributed uniformly in redshift from 
z—0 to 6 (Jakobsson et al. 2006, ). In the simulation we 
assume each GRB releases a radio pulse of luminosity 
1.1 X 10^^ erg s~^ Hz~^ (the luminosity derived in sec- 



tion 4.1.6). In practice, any emission mechanism is likely 
to have a more complex luminosity function, due to a 
range of factors such as initial stellar mass, environment 
and so on. As we have no model luminosity function, we 
fall back to the simplest case of a fixed luminosity. This 
luminosity function also reproduces our overall detection 
statistics. We assume that our detection threshold corre- 
sponds to a flux density of 7 Jy, and that our pulses were 
detected in the range 7 to 10 Jy. This range is somewhat 
larger than the actual range, but nonetheless represents 
a conservative interval. 

The simulation implies we would detect pulses above 
our 7 Jy threshold in 26% of cases (which compares well 
with our observed 25%), but the probability of the flux 
density of a single pulse falling in the observed 7-10 Jy 
range is only 3%. The simulation also implies that a sin- 
gle pulse would have flux density > 10 Jy in 23% of cases. 
Furthermore, by applying the binomial theorem, the sim- 
ulation implies we would obtain our observed statistics 
(2 pulses out of 8 GRBs in the range 7 to 10 Jy) with 



probability 2%. 

This simple simulation suggests that, if the pulses are 
of GRB origin, then the probability of detecting two such 
pulses so close to our detection threshold is small (2%). 
We can say, therefore, that we have a 2% confidence in 
a GRB origin for these pulses. 

4.1.8. Upper limits 

There remains a reasonable probability that neither 
of our pulses are associated with an astronomical phe- 
nomenon. In spite of our best efforts to prove otherwise, 
there could still be some factor we have not considered in 
our telescope configuration that could cause short, dis- 
persed pulse in our data. The fact that the pulses are 
detected at relatively low significance (the brightest is at 
6.6cr), and that no theory has been proposed to explain 
the brightness or time delay of the pulse, implies that 
further observational and theoretical work is desirable to 
confirm such pulses as having a GRB origin. 

We can say with certainty that we did not detect any 
pulses at a significance > 7a. Given our system equiv- 
alent flux density and processing method, we obtain an 
upper limit on the detection of single dispersed radio 
pulses at 1.4 GHz between 200 and 1800s after a GRB 
of 1.27W-1/2 jy^ ^-^Qj.^ g_4 X 10-5 s < u; < 32 X 10"^ s is 
the pulse width. 

4.1.9. Comparison with previous surveys 

Assuming our detections are real, our 7 Jy pulse with a 
spectral index of zero is consistent with the non-detection 
in the all-sky survey of Katz et al. (20031 at 611 MHz, 
and also consistent with the non-detection of the foUow- 



up ob servations of Koranyi et al. ( 1995 ) Dessenne et al 
mM at 151 MHz. 

Our data does not constrain the spectral index well, 
and if we are observing the same mechanism as the 
LB, the spectral index could be as steep as the LB 
{a — —4 ±1). If we assume the LB spectral index, then 
our 7 Jy pulse is still consistent with a non-detection in 
the [Katz et aL| ( 2003 ) survey, but if we scale our result 
to the parameters ot Dessenne et al. ( 1996| ) (the most 
sensitive of the 151 MHZ observa tions) our pulse would 
have a flux density of 57 kJy. The Dessenne et al. ( 1996 1 
observations do not exactly match ours, however. While 
the bandwidth smearing due to dispersion is negligible 
(^ 1 ms for a bandwidth of 700 kHz), the lower time 
resolution of 1.5 s reduces their sensitivity to our 25 ms 
pulse from 16 Jy to 60 Jy. Nonetheless, our pulse should 
easily have been detected by thei r experiment, had they 
observed it. The non-detection of Dessenne et al. (19961 
can be explained if only a fraction ot GRBs emit radio 
pulses. Alternatively, given that the LB spectral index is 
not very well constrained, a modification of the spectral 
index to a = — 2.3 is sufficient to put our pulse at the 



Dessenne et al. ( 1996 1 detection threshold, after account- 
ing tor the lo wer time resolution. Th erefore, it is possible 
to explain the Dessenne et al. ( 1996 1 non-detections if ra- 
dio pulses from GRBs have a somewhat flatter spectrum 
than measured for the LB. 

Conserv atively comp a ring o ur limits with the all-sky 
result by Katz et al.| (20031 of 27 kJy, and assum- 
ing a = —4 as for tne LB, then our limit scales to 
29.4i(;-i/^ Jy at 611 MHz, which is sti ll substantially b et- 



ter. Comparing to the results ofjDessenne et al. ( 1996 ) at 
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151 MHz, our result limit scales to 1.3 x 10^ w Jy as- 
suming the worst case a = — 4 spectral index of the LB 



The sensitivity of the Dessenne et al. ( 1996 1 survey to 
pulses less than their time resolution of 1.5 s is 2Aw~^ Jy 
for w < 1.5 s (ignoring scattering), so their result is more 
constraining for all time resolutions we probed. 

4.2. Long time integrations 

Ultimately, our experiment was not designed to de- 
tect flux changes over long time intervals. Our feed was 
not equipped with a Dicke radiometer, which would have 
ena bled us to calibrate the gain fluctuations (Section 
3.4 1 over the 1 s integration time, but would also have 
contaminated the high-time resolution studies. In addi- 
tion, we did not regularly observe flux calibrators during 
our observations, as this would have ruined the period- 
icity searches (which require continuous data streams). 
Our detection limits are therefore completely limited by 
the temperature-related gain variations in our system, 
of approximately 5%hr^^. we would consider change of 
20%hr~^ as a candidate, which implies that a change in 
source flux density of 760 Jy would have constituted a 
detection. 

We did not detect any variation that was not immedi- 
ately attributable to interference, either from satellites or 
the Sun. We set an upper limit of a change of 760 Jy on 
any long-duration emission (> 1 s) from a GRBs within 
200-1800 s, at 1.4 GHz. 



We have detected no convincing repeating dispersed 
candidates. We also detect no candidates on timescales 
> 1 s, but our experiment was not primarily designed for 
such detections. Nonetheless, we set an upper limit of a 
change of 760 Jy on any long-duration emission (> 1 s) 
between 200 to 1800 s from our GRB triggers. 

The detection of single dispersed pulses in this experi- 
ment is intriguing. Clearly the next step is to determine 
whether these pulses are related to their GRBs, for which 
the key problem is ruling out RFI, statistical fluctuations 
and other equipment-related sources as the origin of these 
pulses. The simplest future experiment to rule out these 
origins is to use a coincidence detection, by employing 
the same telescope setup at two widely separated sites. 
The fact that some 25% of GRBs may be accompanied 
by a radio pulse, even with our relatively poor sensitivity, 
suggests that sensitivity is not the key factor in this ex- 
periment. Therefore, similar dishes, feeds and backends 
can be used. More important parameters of this experi- 
ment are the short on-source time (preferably < 200 s), 
and a wide separation between antennas. A simultane- 
ous detection of a single pulse at two widely separated 
sites, even at 6(t, would almost certainly rule out RFI 
and statistical fluctuations, and render atmospheric ef- 
fects a very remote possibility. If such a detection were 
to be made, the future would be wide open to probe 
the astrophysical and cosmological implications of these 
phenomena. 



5. CONCLUSION 

We have searched for prompt radio emission from 
gamma-ray bursts at 1.4 GHz, using a robotic telescope 
and a pulsar backend. Our telescope was typically on 
source within 200 s of the gamma-ray trigger. 

We detected single dispersed pulses following two 
GRBs at significances > 6a. Simple statistical argu- 
ments, and a null trial based on randomizing channels on 
existing data, rule out random fluctuations as the origin 
of these pulses at 95% and ~ 97%, respectively. The ar- 
rival times of both pulses also coincide with breaks in the 
X-ray light curves. While high DM and absence of adja- 
cent RFI lend credence to an astronomical origin, weak 
impulsive RFI and atmospheric origins of these pulses 
remain a distinct possibility. A simple population argu- 
ments suggests a GRB origin for these pulses of only 2%. 
If the radio pulses are associated with the corresponding 
GRBs, they could be related to changes in the central en- 
gine, in particular the delayed formation of a black hole 
due to spin-down of a rotationally-supported magnetar. 
The non-detection of radi o pulses by pr evious s urveys 
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of ([Koranyi et aL|r 995 ; Dess enne et al.| [T996 Katz et al 
20051 r can be explained by insutticient sensitivity of thos( 



those 

surveys, a somewhat flatter spectrum than measured for 
the Lorimer burst, or the possibility that only a fraction 
of GRBs emit radio pulses. 

If the single pulse is not related to the GRB, we set an 
upper limit on the flux density of radio pulses between 
200 to 1800 s after GRB trigger of 1.27w"i/2 Jy, where 
6.4 X 10-5 s < w < 32 X 10-3 s is the pulse width. This 
li mit is substantially better than the all-sky 27 kJy limit 



of Katz et al. 



stram mg as t 
(11996). 



( [2003 ) at 611 MHz, altho ugh not as con- 
le limits on two GRBs by [Dessenne et al. 
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